Squeezing superfluid from a stone: 
Coupling superfluidity and elasticity in a supersolid 
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Starting from the assumption that the normal solid to supersolid (NS-SS) phase transition is 
continuous, we develop a phenomenological Landau theory of the transition in which superfluidity 
is coupled to the elasticity of the crystalline ^He lattice. We find that the elasticity does not affect the 
universal properties of the superfluid transition, so that in an unstressed crystal the well-known A- 
anomaly in the heat capacity of the superfluid transition should also appear at the NS-SS transition. 
We also find that the onset of supersolidity leads to anomalies in the elastic moduli and thermal 
expansion coefficients near the transition and, conversely, that inhomogeneous lattice strains can 
induce local variations of the superfluid transition temperature, leading to a broadened transition. 
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Superfluidity — the ability of liquid "^He, when cooled 
below 2.176 K, to flow without resistance [1, 2] 
through narrow pores — has long served as a paradigm 
for the phenomenon of "off-diagonal long-range order" 
(ODLRO) in quantum liquids and superconductors [3]. 
Supersolidity — the coexistence of ODLRO with the crys- 
talline order of a solid — was proposed theoretically [4- 
12] as an even more exotic phase of solid "^He, but it 
has eluded detection [13, 14]. Recently, Kim and Chan 
[15, 16] have reported the onset of "nonclassical rota- 
tional inertia" [6] in a torsional oscillator experiment with 
solid "^He, and they interpret their results as indicating 
the onset of supersolidity. However, their interpretation 
remains controversial [17-22], so it is important to com- 
plement the nonequilibrium torsional oscillator measure- 
ments with equilibrium thermodynamic measurements, 
e.g., of the specific heat. In this work we start from 
the assumption that normal solid to supersolid (NS-SS) 
phase transition is continuous, and develop a phenomeno- 
logical Landau theory of the transition in which super- 
fluidity is coupled to the elasticity of the crystalline ''He 
lattice. Wc find that the elasticity does not affect the 
universal properties of the superfluid transition, so that 
in an unstressed crystal the well-known A-anomaly in the 
heat capacity of the superfluid transition should also ap- 
pear at the NS-SS transition. We also find that the onset 
of supersolidity leads to anomalies in the elastic moduli 
and thermal expansion coefficients near the transition; 
conversely, inhomogeneous strains in the lattice can in- 
duce local variations of the superfluid transition temper- 
ature, leading to a broadened transition. As our theory is 
rooted in a few simple assumptions and symmetry princi- 
ples, we expect our results to be robust and insensitive to 
the details of a microscopic model for the supersolidity. 

We hypothesize that, as with the superfluid ^He, the 
appropriate order parameter describing the onset of su- 
persolidity is a complex scalar fleld V'(r), depending on 
the location r. We shall assume that the phase transition 
to the supersolid state is continuous, as is the superfluid 



"^He transition. Then, as is well known from the theory of 
critical phenomena [23] , the universal properties of the 
supersolid transition may be obtained via a model that 
retains only those terms in the free energy up to leading 
(relevant) order in powers of the fields and their spatial 
gradients, resulting in the Landau form: 
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(1) 



Here, a{T) depends smoothly on the temperature T; it 
is negative at low temperatures, changing sign slightly 
above the transition temperature T^. [If fluctuations of 
"0 are ignored, would be exactly the temperature at 
which a(T) changes sign.] For T < (T > T^) the 
thermal expectation value of (ip) is non-zero (zero). The 
constant w measures the strength of the nonlinearity, 
and the symmetric tensor Ca0 characterizes the spatial 
anisotropy inherited from the crystallinity of the normal 
solid. For an isotropic superfluid or for a cubic crystal 
Cq/3 = cSap, with c a constant, whereas for an hep crys- 
tal (such as solid helium [24]) is uniaxial, such that 
Cq/3 = CzTLanp -\- c_l((5q/3 — nanp), with n a unit vector 
that points along the preferred axis of the crystal, and Cz 
and c_L independent constants. It is important to note 
that the symmetry of the superfluid density tensor 
(which relates the superfluid velocity to the momentum 
density) is the same as that of c^^. 

If we were dealing with the normal-to-super/?wirf tran- 
sition, Eq. (1) (with Cap = cSafs) would be the entire 
story. However, because we are dealing with the normal- 
to supersolid transition, it is not: when formulating the- 
ories of continuous phase transitions, it is necessary to 
keep all of the degrees of freedom that are soft (i.e. ex- 
hibit large thermal fluctuations) at the transition [23]. 
For a superfluid, is the only such degree of freedom. 
However, in a supersolid there are additional phonon de- 
grees of freedom that are soft not only at the transition 
but throughout both the normal and supersolid phases. 
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These are associated with displacements u(r) of the po- 
sitions r of the imdistorted crystal lattice. (The normal- 
solidification of the ^Hc occurs at a temperature substan- 
tially higher than so that amplitude fluctuations of the 
density waves are not soft and may be neglected.) 

As the free energy must be invariant under spa- 
tial translations and rotations [25], it can only de- 
pend upon u(r) through the symmetric strain tensor 
Uai3 ^ ^ {daU/3 + d/sUa + daU^ djsUj). Thus, at long 
wavelengths the relevant terms in the free energy involv- 
ing Ua/3 alone are simply those of standard elastic the- 
ory: J-c = ^ J d^r Xap-ySUapUys, where Xap-iS arc the 
bare elastic constants, and repeated indices are summed 
over. The form of Xaff-yS is dictated by the symmetry of 
the crystal; for an hep crystal such as ''He, it is parame- 
terized by 5 independent elastic constants [25]. 

To determine the form of the coupling between the su- 
persolid order parameter ip and the local displacement 
field u(r), we follow the work of Aronovitz et al. [26] 
by allowing the (formerly constant) parameters in J^g 
[i.e. a{T) and w] to depend on the local value of u in 
a manner consistent with the symmetries of the system. 
Thus, we expand in powers of the strain tensor, a{T) — > 

a^°^ + a-^l u^p + a[fj^,_5 u^^ u^g H , and similarly for 

w. Here, the tensor aa/s has the same symmetry as Ca/3, 
and hence for a uniaxial crystal (such as the hep phase of 

solid helium [24]) a^^^ = a^Uarip + a±{Saf3 — ^^c^?^/3)^ with 
Gz and a_L independent coupling constants. 

The terms of 0(u^^) in the expansion of a and of 
0{uap) in the expansion of u; all prove, by nai've power 
counting [23], to be irrelevant; i.e., they do not affect 
the universal critical properties of the transition. Hence, 
these properties follow from the following minimal model: 

+ 2Aa/3^i W75 + 2««i"a/3 1 V^r} • (2) 

In our minimal model the effect of the elasticity is to pro- 
duce a local, strain- dependent critical temperature for the 
superfluid. In fact, our minimal model is formally equiv- 
alent to that of a magnetic system of planar spins [an 
0(2) ferromagnet] on a compressible lattice, for which 
locally dilating or compressing the lattice causes the ex- 
change couplings to decrease or increase, resulting in a 
local change of the critical temperature. With this anal- 
ogy in hand, we may use the work of de Moura et al. [27], 
which shows that the elasticity is irrelevant to the critical 
properties of the 0{2) ferromagnet, provided the specific- 
heat exponent a of the decoupled (i.e. a^^'j = 0) system 
is negative, which it is at d = 3 [28]. Thus, we may 
conclude that the universality class, and hence univer- 
sal properties, associated with the supersolid transition 
are unaltered by the coupling of elastic degrees of free- 
dom to the supersolid order parameter. In particular, 
this implies a A-anomaly in the specific heat Cp near the 



transition: 

CpW = — |i|-"(l + + 6^ + •••)+ (3) 

a 

where t = {T — Tc)/Tc denotes the reduced temperature, 
a = — 0.0127 ± 0.0003 is the universal specific heat expo- 
nent of the superfluid transition, and A = 0.529±0.009 is 
the equally universal correction to scaling exponent [28] . 
The subscripts -|- and — in Eq. (3) denote behaviors above 
(T > Tc) and below (T < Tc) the transition, respectively. 
Although the constants A± are non-universal (they will, 
e.g., change as one moves along the SS-NS phase bound- 
ary in the pressure- temperature phase diagram), their 
ratio is universal; the current best estimate of its value 
is = 1.053 ±0.002 [28]. The superfluid density 

tensor also exhibits universal scaling with reduced 
temperature t, i.e., = °/o^^ where the tensor °p^^ 
is temperature independent, and jy = 0.67155 ± 0.00027 
is [28] the universal correlation- length exponent. 

Despite the irrelevance of the elastic couplings to ip for 
the universality class of the transition, these couplings do 
have important, experimentally observable consequences: 
due to them, the elastic properties inherit singularities in 
their temperature dependence from parent singularities 
associated with the critical fluctuations of the supersolid 
order parameter. Following Ref. [26], we construct the ef- 
fective free energy J^^s governing the elastic fluctuations 
by integrating out the supersolid fluctuations: 
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Proceeding perturbatively, in powers of the tp-Uap cou- 
pling term, we obtain (neglecting an additive constant) 



(5) 

where the singular parts of the fluctuation corrections to 
the bare elastic constant and stress tensors are given, to 
leading order, by 



(6a) 
(6b) 



Here, the governing functions, C (which proves to be 
proportional to the singular part of the specific heat) and 
£), are given by 
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C(t) = lim 



k^o J (27r)'* \ 2 2 



(7a) 
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where the expectation values (• • • )o are taken with re- 
spect to the rigid ■0 system (i.e. a system in which all 
Uafj's are frozen at zero), whose free energy ^rig is just 



3 



J^csiua/s] with all u^p's set to zero. Moreover, \^p\'^ de- 
notes the Fourier transform of |'i/'(r)p, and (• • • indi- 
cates a connected correlator. In addition, j^^jg is the free 
energy density associated with J^ng- 

D{t) and C{t) can readily be related to the specific 
heat singularity of the rigid system, by noting that a^^^ 
is a linear function of T close to T^, so that derivatives 
with respect to a'-"-' are proportional (near T^) to deriva- 
tives with respect to T. Hence, simple thermodynamic 
identities imply that D(t) and C{t) arc proportional to 
the entropy and the specific heat of the rigid system, re- 
spectively, near Tc- This immediately determines their 
critical behavior: D{t) = G±|i|i-"(1 + a±|t|'^ + ■■■) 
and C{t) = -ha±|t|^ + ■■■), where G+/G" = 

A'+/A'~- = A+/A- ^ 1.053 ± 0.002 are universal. 

It is evidently valuable to estimate the size of the 
expected elastic and thermal expansion anomalies. To 
do this, we need an estimate of the couplings (^^^ in 
Eq. (2), and this can be obtained, following Ref. [29], 
from the form of the SS-NS phase boimdary, Tc{P), in 
the pressure-temperature phase diagram. In particular, 
the couplings are related to the slope of the boundary via 

aTe(P)/aP~a«/(Aa'), (8) 

where a' = da'-"'^ /dT, and A is a typical elastic constant. 
With this in hand, we can now estimate the size of the 
elastic and thermal expansion anomalies. This can be 
accomplished by using Eq. (7) to estimate C and D well 
away from the critical point (say, at T = 2Tc), where we 
can make the Gaussian approximation to the correlation 
functions, which gives, e.g., 
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where is the high-temperature correlation length for ij) 
fluctuations. We have taken the integral over k to have 
an ultraviolet cutoff comparable to ^o, used the fact that, 
well above Tc, a^^\T)/c ~ ^q"^ to replace the propagator 
[up to 0(1) factors] by a^^^ for all fc, and estimated a{T = 
2Tc) ~ a'Tc- For supeyfluid ^He, the length is known 
to be comparable to the atomic size: ~ 0.2 nm; for 
want of better information, we shall assume that this is 
also true for the supersolid. 

By using this estimate for D{T) in Eq. (6b), and then 
minimizing Eq. (5) over u„/3j we arrive at a typical value 
for the thermal expansion: 



Su, 



■a/3 



d fcBTc(P) 
dP ^1 ■ 



(10) 



where in the last step we have used Eq. (8). Estimating 

\dTc{P)/dP\ - 10-2 K/atm by its value in the liquid 
state of ''He [24] gives a typical value of Su^js ~ 0.16. 
For the fractional anomaly in the elastic constants, argu- 
ments essentially identical to those just used give 
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0.17. 



(11) 



Perhaps the best way to observe the predicted anoma- 
lies in the elastic constants is through sound speed mea- 
surements. In single crystals, the polarization of the 
sound modes studied must be chosen judiciously: the uni- 
axial form for the expansion coefficients, a^^^ = aznani3 + 
«-L(^a/3 — naUf)), implies that only "bulk" elastic mod- 
uli, specifically, the terms ^zzinanpUafiY' + Aj^j^(((5a^ — 
nani3)uapY + A_Lz(((5a/3 - nani3)uap){nanj3Uaf3), acquire 
anomalous temperature dependence. The elastic con- 
stants A22, A_L_L, and are readily shown by a stan- 
dard sound mode analysis [30] to only affect the sound 
speeds of modes polarized in the plane formed by the 
normal to the hexagonal layers and the direction of prop- 
agation. Furthermore, in this plane, transverse modes 
propagating cither along, or orthogonal to, the layers 
also have sound speeds independent of the anomaly- 
displaying elastic constants \zz, Ax±, and \_\_z. Hence, 
to observe the anomaly in the sound speeds in a single 
crystal of supersolid hep helium, one should study modes 
polarized in the plane formed by the hexagonal layers and 
the direction of propagation, and choose that propaga- 
tion direction not to lie in, or orthogonal to, the layers. 

It is, however, unlikely that experiments will be per- 
formed on single crystals of helium. It is far more likely 
that they will be performed on polycrystalline samples, 
which are macroscopically isotropic, due to the random 
orientations of the constituent crystallites. Calculating 
the isotropic shear modulus ji and Lame coefficient A of 
such an ensemble of randomly oriented crystallites is well- 
known to be a formidable problem (see, e.g., Ref. [31]). 
Nonetheless, using the exact bounds of Hill [32], we can 
show [30] that both the shear modulus and the bulk mod- 
ulus of a macroscopically isotropic polycrystalline helium 
sample will exhibit the \t\~°' anomaly we predict here. As 
a result, both the transverse and the longitudinal sound 
speeds of such a sample will show the [t]^" anomaly. 

Beyond the critical properties described above, our 
model has the important implication that in a he- 
lium crystal the SS-NS transition would be rounded by 
any spatially inhomogeneous internal stresses that make 
'"a/3 7^ 0. As such stresscs are almost unavoidable in 
any crystal (and are believed to be present in the exper- 
iments of Kim and Chan [15, 16]), this roimding is al- 
most certain to be present in all experiments performed 
to date. The reason for such broadening is very sim- 
ple: stresses that make a|^^^ spatially inhomogeneous 
make of the SS-NS transition spatially inhomogeneous 
as well; cf. Eq. (2). Therefore, roughly speaking, distinct 
parts of the sample would become supersolid at distinct 
temperatures, broadening the transition. Such broad- 
ening is evident in the p^{T) vs. T plots of Kim and 
Chan [16], which clearly do not show the expected \t\'' 
singularity near the putative Tc, with v being the corre- 
lation length exponent. [If they did, p^{T) vs. T would 
hit the horizontal axis perpendicularly, rather than — as 
in those data — tangcntially.] This broadening may also 
explain the apparent absence of the expected A-anomaly 
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in the specific heat [33] near the putative Tc in those 
experiments: this peak is simply "smeared away'.' 

Our model suggests that the best way to definitively 
observe (or rule out) the critical behavior of the specific 
heat, elastic constants, superfluid density, etc. that we 
predict here, and thereby to test the notion that helium 
actually does exhibit a supersolid phase, is to prepare the 
solid helium samples in a way that is designed very care- 
fully to eliminate any spatially inhomogeneous stresses 
in the crystal. At the very least, this would require that 
samples always be pressurized at high temperatures in 
the liquid phase, and only then cooled into the solid 
phase; this cooling should be done at constant pressure, 
so as to avoid introducing inhomogeneous strains. 

One intriguing, albeit highly speculative, final impli- 
cation of the coupling of strains to the supersolid order 
parameter has to do with the very existence of the 
supersolid state in "^He. A number of microscopic cal- 
culations [17] suggest that ^He does not have a super- 
solid state at all, in contradiction with the experiments 
of Chan et al. [15]. These calculations were presum- 
ably done under conditions of purely hydrostatic stress 
(i.e. simple pressure), for which only the three diagonal 
components of the strain tensor are non-zero and equal 
(i.e. Uxx — Uyy — Uzz 7^ ) . If it were the case that the 
coefRcients and a± in the uniaxial coupling tensor for 
solid hep '*He happened to obey Uz « —2a± then the 
effect of this pressure on the effective Tc for supersolid 
order would be very small. If the unstrained Tc were 
negative [i.e. the coefficient a^°' in Eq. (2) were positive 
for all T], this would imply that, under such a hydro- 
static stress, the crystal would never enter the super- 
solid state. If, however, the crystal were to be subjected 



to an anisotropic stress (i.e. one for which the relation 
Uxx = Uyy = Uzz is not satisfied), the near cancellation of 
a^^j "a/3 would not occur, and this term might be able to 
make the effective Tc positive [i.e. change the overall sign 
of the coefficient of in Eq. (2)]. That is, it is possi- 
ble in our model that, although a hydro statically stressed 
sample would not show a supersolid phase, an anisotropi- 
cally stressed one might. This could be highly significant, 
as Chan et al. [15] believe that such anisotropic stresses 
are present in their samples. Although the above ar- 
gument is obviously quite speculative, it seems to us, 
nonetheless, a possibility that these random stresses, far 
from being an experimental nuisance, might just be what 
is responsible for the presence of supersolidity. Support 
for this idea comes from two other facts: (1) Many ex- 
periments [13, 14] do not see supersolidity. Perhaps these 
samples simply lacked sufficiently large inhomogeneous, 
anisotropic stresses. (2) The superfluid fraction in the 
experiments of Chan et al. [15] is extremely low. Could 
this be because only very small, highly anisotropically 
stressed regions of the sample are going supersolid? 

In summary, by assuming that the NS-SS transition is 
continuous and using general symmetry arguments, we 
have constructed a minimal model for understanding the 
NS-SS transition. Our model predicts a rich phenomenol- 
ogy, whose existence, if observed, would help corroborate 
the torsional oscillator evidence for the supersolid phase. 
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